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The super-inflationary phase is predicted by the Loop Quantum Cosmology. In this paper we 
study the creation of gravitational waves during this phase. We consider the inverse volume correc- 
tions to the equation for the tensor modes and calculate the spectrum of the produced gravitons. 
' The amplitude of the obtained spectrum as well as maximal energy of gravitons strongly depend on 

\ the evolution of the Universe after the super-inflation. We show that a further standard inflationary 

. phase is necessary to lower the amount of gravitons below the present bound. In case of the lack of 

OA ' the standard inflationary phase, the present intensity of gravitons would be extremely large. These 

considerations give us another motivation to introduce the standard phase of inflation. 
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I. INTRODUCTION 



The cosmological creation of the gravitational waves was proposed by Grishchuk in the mid-seventies. Since 
that time this phenomenon has been studied extensively, especially in the context of the inflation. The accelerating 
expansion phase gives the conditions for the abundant creation of the gravitational waves. Gravitons produced during 
the inflation fill the entire space in the form of a stochastic background. Together with the scalar modes, produced 
during the inflation, they form primordial perturbations leading to the structure formation. The analysis of the cosmic 
~^ microwave background (CMB) and large scale structures gives therefore the possibility of testing inflationary models. 
^ In the case of the CMB the impact of the gravitational waves comes from the primordial spectrum and from tensor 
^ I Sachs- Wolfe effects. The Sachs- Wolfe effect is somehow secondary and leads to the CBM anisotropics as the result of 
the scattering of CMB photons on the relic gravitons. The form of this anisotropics is given by 
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, where /ly describes tensor modes and is the vector parallel to the unperturbed geodesies. The influence of the 
gravitational waves for the CMB is however to weak to be observed directly with the present observational abilities. 
Another possible method to detect gravitational waves is to use of the antennas like LIGO, VIRGO, TAMA or GEO600 
■ 01 @| • Although these detectors are now very sensitive this is still not enough to detect directly the gravitational 
$H ' waves background It may look pessimistic, we hope however that some further improvement of the observational 
skills bring us the observational evidence, so needful for the further theoretical improvements. 

In this paper we consider a new type of the inflation which naturally occurs in the Loop Quantum Cosmology 
. This is so called the super- inflationary scenario 0, Q and is a result of the quantum nature of spacetime in the 
Planck scales. The spacetime is namely discrete in the quantum regime and its evolution is governed by discrete 
equations. However for the scales greater than a; = y^lp\ (7 is so called Barbero-Immirzi parameter) the evolution 
of the spacetime can be described by the Einstein equations with quantum corrections. For typical values of the 
quantum numbers, super-inflationary phase takes place in this semi-classical region. 

Our goal is to describe the production of the gravitational waves during the super-inflation. This problem was 
preliminary analysed in Ref. [8|, but quantum corrections to the equation for tensor modes was not included to 
calculate the spectrum of the gravitons. In this paper we include so called inverse volume corrections to the equation 
for evolution of the tensor modes and then calculate the spectrum of produced gravitational waves. The equations 
for the tensor modes was recently derived by Bojowald and Hossain They had analysed the inverse volume 
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corrections and corrections from holonomies. In this paper we concentrate on these first ones. The quantum corrections 
are generaUy comphcated functions but they have simple asymptotic behaviours. To calculate the productions of 
gravitons during some process we need somehow to know only initial and final states, where asymptotic solutions are 
good approximation. In these regimes calculations can be done analytically. We use numerical solutions to match 
them. 

The organization of the text is the following. In section II we fix the semi-classical dynamics. Then in section III 
we consider creation of the gravitons on the defined background. In section IV we summarize the results. 



II. BACKGROUND DYNAMICS 

The formulation of Loop Quantum Gravity bases on the Ashtekar variables and holonomies. The Ashtekar 
variables replace the spatial metric field qab in the canonical formulation as follow 

a: = K + ^Kl (2) 
= (3) 

where FJ^ is the spin connection defined as 

n = -e^^-^-4(9[,e^] + le^^eiVi]) (4) 

and the Kl is the intrinsic curvature. The is the inverse of the co-triad ejj defined as Qab = e^e^. In terms of the 
Ashtekar variables the full Hamiltonian for general relativity is a sum of constraints 

i/*?' = J d^x (iV'G, + N'^Ca + NK,), (5) 



where 



Ca = E^Fl, - (1 ~ i^)KG,, 

G, = DaEt (6) 



and the scalar constraint has a form 



Hg ■■= J d^xN{x)hs, 



' d'xNix) ( ^,I=2=e^hF^, 2{1+ j')^=^KlKi (7) 



16^G J ' ' y VT^etE] ' ' ' ^\AetE\ 

with F — dA+ i [A, A\ . The full Hamiltonian of theory is a sum of the gravitational and matter part. With convenience 
as a matter part we choose the scalar field with the Hamiltonian 

We assume here that field (j) is homogeneous and start his evolution from the minimum of potential V{(j)). The second 
assumption states that contribution from potential term is initially negligible. So the density of Hamiltonian is 
simplified to the form Ti^ = (l/2)7r^/-\/| det E\. The term 1/ det E\ for the classical FRW universe corresponds to 
where a is the scale factor. On the quantum level term deti?| is quantised and have discrete spectrum. In 
the regime a ^ we can however use the approximation det E\ = D /a^ where 

2 2/) 

D{q) = g^/^ {I [{q + -\q- 1^'] - ^^ii^i + 1)'+' - ^^iq - l)k - } ' 0) 

and q = (a/a^,)^ with a, = a/tT/S^pi- Function ([9]) depends on the ambiguity parameter I. As it was shown by 
Bojowald [ll| the value of this parameter is quantised according to = 1 — (2fc)^^ > 1/2, fc € N. For the further 
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investigations we choose the representative value Z = 3/4. In the semi-classical region a, ^ a ;:!> a.; expression ^ 
simplify to the form 



where 



D = D^a" 



o X 3/(2-2/) 

^ a^^^^"'^/^^"'^ and n 3(2 - 0/(1 - 
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(10) 



(11) 



Now, due to the Hamilton equations we can derive the Friedmann and Raychaudhuri equations for the flat FRW 
universe filled with a homogeneous scalar field 



H' = 

a 

a 
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The equation of motion for the scalar field with quantum corrections has the form 



(12) 
(13) 

(14) 



As we mentioned before, for the further investigations we simplify equations p2|) . (fT3l) and (|14p assuming V{(j)) = 0. 

The expression for the quantum correction D is complicated and it is impossible to find an analytical solution for 
the equations of motion. In fact we even do not need it for the future investigations. To calculate the spectrum of 
gravitons we need to know analytical solutions only for the inner and outer states. We choose the |in) and |out) states 
respectively in the quantum and classical regimes. The expression for the quantum correction ([9]) simplifies to the 
form D = D^a^ for the a; < a <C a* and D = 1 for a ^ a*. In these limits we can find the analytical solutions for the 
equations of motion (fT2|) . (fT3|) and (fT4| . It is useful to introduce the conformal time dr = dt/a to solve equations and 
for the further investigations. In the next step we must to fit obtained asymptotic solutions using a global numerical 
solution. The solution for the evolution of the scale factor in the quantum limit has the form 




where p = 2/(4 — n). The solution in the classical limit we obtain putting simply I = 2 what gives D 
The constants of integration ^ and f3 we fix with the use of a numerical solution applying formula 



P a| 



and P = — T — p—f- 

n' 



(15) 

1 andp = 1/2. 



(16) 



The value of the conformal time r must be chosen in a proper way for the given regions. We will discuss this question 
in more details later. 

Our point of reference is the numerical solution. To make this description complete we must choose the proper 
boundary conditions for the numerical solution. We use here the condition for the Hubble radius which must be larger 
than the limiting value ai [l2|, what gives us 
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\H\a<- 

a,. 



\H\a^ < 1. 



(17) 



The next condition requires that the scale factor must be greater than at the bounce. It is fulfilled taking a\rg = a* 
for some value of the conformal time tq. In fact the conformal time is unphysical variable and their value can be 
chosen arbitrary. The physical outcomes do not depend on coordinates because the theory is invariant under local 
diffeomorphisms. So as an example we can choose 



=-4 = a* 



\to = -4 



The chosen value of a'\Tg^_4 holds the condition ([T7]). Namely for j = 100 we have \i 
solution is shown in Fig. [T] as a black line. We fix the boundary approximations in ri 



: 0.084 < 1. 

-20 and T2 



(18) 
(19) 

The numerical 
= -1. The the 
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FIG. 1: Numerically calculated evolution of the scale factor (top black line) and the Hubble parameter (bottom black line) 
for the model with j — 100, / = 3/4 and with initial conditions (|18|l and (|19p . The approximated solutions described by psp 
and ((22)) . 



numerical solution gives us for these points 

a\r,=-20 = 0.536?Pi , a'|^,=_20 = 0.007?Pi (20) 
a\r,_=-i ^ 3.912?pi , a'|^,=_i = 0.461/pi. (21) 

Now with the use of expressions (jl6p we can fix the approximated solutions. However when we use formula (jl6p 
directly to calculate parameters in solution (|15p for the outer state we obtain complex ^. It is due to the expression 
under square {p ~ 1/2) is negative. So to put away complex numbers we redefine the exit solution to the form 

a = K^T + C (22) 

where 



K = -ze (23) 

C = -/?• (24) 

We show in Fig. [1] how these approximated solutions match with solutions obtained numerically. As we can see 
these approximated solutions well describe the evolution in the neighbourhood of a*. 



III. GRAVITATIONAL WAVES 



We have already mentioned in section|T]that gravitational waves can be abundantly produced during the accelerating 
phase. In this section we want to show in details how it works and calculate properties of produced gravitons. To 
describe the spectrum of gravitons it is common to use the parameter 

= (25) 

Pc av 

where pgw is the energy density of gravitational waves and pc is present critical energy density. Our goal in this 
section is to calculate the function U,g-„[v) for the gravitons produced during the super- inflationary phase. 
The gravitational waves hij are the perturbations of the background spacetime in the form 

ds^ = a^{T) [-dT^ + (% + h,j)dx'dx^] (26) 

where \hij \ <C 1. Using constraints h\ = Vi/i* = we can see that tensor have only two independent components 
h\ = —h\ = and h\ = h\ = h-^. These components correspond to two different polarisations of gravitational 
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waves. Inserting the perturbed metric ([26)) to the Hilbert-Einstein action S'h-e = (l/167rG) / d'^x^/—gR gives the 
series ^h-e = S^^^ + S^^^ + 5*'^' + . . . , where the second order term have a form 



1 



647rG 
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32nG 



(27) 



and give us the action for the gravitational waves. The two kinds of polarisations are not coupled and can be treated 
separately. To normalise the action and simplify the notation it us useful to introduce the variable 



h = 



Vl67rG Vl67rG' 
what leads to the expression for the action in the form 



(28) 



(29) 



This action is the same like the action for an inhomogeneous scalar field without the potential. Inverse volume 
corrections can be therefore introduced in the same way like in the case of the scalar field. As we mentioned in 
Introduction there are also holonomy corrections to this action. Here we consider however the influence from the 
better examined inverse volume corrections. 

As it was shown by Mulryne and Nunes [l^, in the context of scalar field perturbations, it is useful to introduce 
the variable u = aD~^/^h and rewrite the action (j29p with quantum corrections to the form 



(30) 



where 



'cff - 



D 

a 



(31) 



Till now the considerations of the gravitational waves has been purely classical. The next step is the quantisation of 
the classical gravitational waves what brings us the concept of gravitons. To quantise the field u we need to firstly 
calculate conjugated momenta 



(32) 



The procedure of quantisation is the simple change of fields u and tt for the operators just adding hats and to introduce 
the relations of commutation. We decompose operators considered for the Fourier modes 



u{t, x) 
7r(T, x) 



1 



2(27r)3/2 
1 



MT)e'^-'' + ul{T)e 
7rk(T)e^'^-'' + 7rt(r)e- 



— zk-: 



■ikx 



2(27r)3/2 

where the Fourier components fulfil the relations of commutation 

[uk(r),7rt(r)] = *<S(^)(k-p), 

[uk(T)t,7rp(r)] = tS^^\k~p), 

["k(T),*p(T)] = *<5(3)(k + p), 

["k(r)t,7rt(r)] = zS^^\k + p). 



(33) 
(34) 



(35) 
(36) 
(37) 
(38) 



To express the Fourier modes in terms of the annihilation and creation operators we need to solve the quantum 
Hamilton equations 



u 

TT 



i[Ht,u], 

l[Ht,iT]. 



(39) 
(40) 
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The Hamilton operator have the form 

A. ' 



1 1 



2 4(27r)3 
DS'H 



+ m 



cS 



ik-x 



VTke 



tg-jk-x 



TTkTT^ + TT^VTk + (iJfc^ + m^g) ( Mk^k + "k'^k 



(41) 



where we inserted decompositions ([SSj) and ((34)) . When we apply the Hamiltonian (|4T|) and the decompositions ([33] 
and p4p. the Hamilton equations (|5^ and (|l(I)l take the forms 



-^k = T^k, 

Tk = - + mlg) Uk- 



(42) 
(43) 



The general solution of these equations has the form 

Mr) = ak/(fc,T)+aLkr(fc,T), (44) 
7rk(T) = ak5(A:,r) + a^i^g*(fc,T). (45) 

where /(fc, t)' = g(k, t). When we insert these solutions to the Fourier decompositions (I55|) and (IM)) we simply obtain 

1 



u(t,x) = 
7r(T,x) = 



(27r)3/2 
1 



g{k, T)d]^e'^''' + g*{k, T)aj,e" 



(27r)3/2 

The mode functions fulfils the so called Wronskian condition 

r{k,T)g{k,r)-f{k,T)g*{k,T)^-t 



(46) 
(47) 

(48) 



as a result of relations of commutation (|35ll38p . These relations is important to normalise properly the mode functions. 
The Hamilton equations (|^^ and ([33]) together with give us the equation for the mode function 



(49) 



dr^ 



f{k,T)+[De + mlJfik,T)^0 



This equations has two regimes. The first one called adiabatic corresponds to the situation when Dk^ + m^ff = F ^ 0. 
The second one leads to the super-adiabatic amplification and corresponds to the situation when F <C 0. The creation 
of the gravitational waves corresponds to the case of the super-adiabatic amplification. 

We can now investigate which modes are amplified. The condition F ^ corresponds to Dk^ <C —mlg. In Fig. [2| 
we see the evolution of —rri^g. As we see the condition for the creation of gravitational waves is fulfilled in the region 
of the super-inflation (r G [~ — 6,~ —4]) . In the right panel we can see the evolution of F for three different values 
of fc. 

The value of k is dimensionless in the undertaken scheme. To obtain a dimensional value we must multiply it simply 
by the corresponding scale factor which has a dimension of length. The wave number k, for example for the final state 
at, corresponds to the frequency 



_af_\ 
27raf \atoday J 

measured today. The important task is to calculate the factor atoday/of- We can make the decomposition 

^today O^today ^dec '^rh ^i-ciid *^i-start 



Of 



'^'dcc O-rh Q-i-cnd Q-i-sta 



af 



(50) 



(51) 
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FIG. 2: Left : Evolution of the parameter —rri^g for the model considered. Right : Evolution of the function F — Dk^ + 
for different values of k. 



where 

fldoc —> photons decouphng, Zdcc — 1070, Tdcc — 3 • 10"^ K — 0.2 eV 

Orh — *■ reheating phase , T^h = Tgut — 10^^ GeV 

ii-ond — ^ end of inflation 

fli-start beginning of inflation 

We can assume the sudden reheating approximation (a^h — ai_cnd) and the standard value of e- folding number for 
inflation N = ln(ai_cnd/ai-star) = 63. Then increase of the scale factor forms the final state till present value assumes 

~ (1 + zdoc) • ^ • • ~ 10=^^ (52) 

We had assumed here that ai_start/af — 10^, this value can be obtained from numerical simulations like these in 
Ref. [14]. Now we can return to equation ([50)1 and then calculate the maximal frequency of produced gravitons. 
From the right panel in Fig. [2] we can see that in the vicinity of r ~ —6 we have the transition from positive to 
negative values of F. So as we can see for some 10 < fck < 1000 we have the transition from the adiabatic to the 
super-adiabatic regime. From the numerical investigation we obtain fctr — 500. In fact the function F, for small 
values of r, is always slightly below zero even for greater values of k than fcti-. But we assume that this effect is 
negligible. In fact higher values of k easily reach the transplanckian scales and it is not clear that we should trust 
the standard physics in this regime. The wave number ktr corresponds roughly to the scales a,. So the value of ktr 
corresponds to k{ = (a*/af)fctr — 250 in the final state ai{T = — 1) = 3.912Zpi. Applying equation ([50]) we obtain the 
maximal frequency for the present epoch //max — 2 • 10~^^ Hz. This value corresponds to the scales Amin — 5 kpc. It 
is instructive to consider the model without inflation. In this situation the maximal frequency of the relic gravitons 
would be i^max — 2 • 10^^ Hz what is extremely huge number. In fact it is possible that GUT energy scale and inflation 
cover each other in some place. It was somehow one of the motivation to introduce inflation to solve the problem of 
topological defects. In the case when the GUT scale occurs after thee reheating the problem of topological defects 
must be solved in a different way. We mention this problem to show that it is possible that the value Otoday/if can 
be lower than calculated before. 

After this analysis we can return now to equation (|49p . We need to use the deflnition of the quantum correction 
D and effective mass rn^g in the quantum regime. It is also useful to rescale the conformal time and introduce 
—7] ~ —T + (3. We will back later to the previous definition because an additional degree of freedom (3 is necessary to 
fit properly the boundary solutions with numerical one. Equation (I49p takes the form 



-^fik, r,) + [i?4"(-^)"Pfc2 - p{p - 1)] f{k, r,) = (53) 

and the general solution in terms of Bessel functions have the form 

f{k,ri)^Ci^J\,\[x)+C2^Y\,\{x) (54) 
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with 



\l + np\ 

Vl + 4p(p-l) 
^ = • 

2 + np 

With the use of the Wronskian condition (l48l) we can rewrite the sohition (l54l) to the form 

/(fc,'7) 



2|2 + np 

where we introduced Hankel functions defined as 



-'7 



D,H\;^l{x) + D^H\:;l{x) 



(55) 
(56) 

(57) 

(58) 
(59) 



and the constants Di and D2 enjoy the relation \Di\'^ — I-D2P = 1- To fix values of the constants D\ and D2 we must 
consider the high energy limit, namely a; ^ 1. fn this limit the Bessel functions behave as follow 



J\v\{x) ■ 



sm X 



nx 
2" 



~2 4 

|t^|7r TT 

~2 4 



(60) 
(61) 



what give us H^^^^{x) ^/2/{^^x) exp [ix — i\v\'K /2 — iir /4\ and H^^'l{x) ^2/{'kx) exp [—ix + i\v\'K /2 + iir Clas- 
sically the limit cc ^ 1 corresponds to advanced solution called the Bunch-Davies vacuum e~*'^'^/-\/2fc. Generally the 
limit obtained here differs from the classical one but can be restored taking 1 — 2. Then to obtain the proper high 
energy limit we must choose D2 = and Di = exp [i|i/|7r/2 + in/A]. Applying evaluated values of Di and D2 to the 
solution ()57p we finally obtain mode functions for the initial state 



Mk,r) = Af^^^-kT + kf3H\l]{ 



k 



-kr + k(3 



1 h\1\{x) 2 + 



np 



2 -kr + kl3 



■np\ 



where 



e'^(^+f) 



2\2 + np\ 



(62) 



(63) 



(64) 



The functions gi{k, r) were calculated from relation f{k, t)' = g{k, r) showed earlier. We had used also the expression 
for derivative of the Hankel functions in the form 



(65) 



Similar investigations lead to the expression for the mode functions for the final state. We must remember that it is 
however not only a simple set ^ = 2 but also the change of the solution for the scale factor to this expressed by \22\ . 
In this case we obtain 



IfiKr) = e^^ JJL ^fk7TKH^^\kT + kO, 



I Trk n — 

— v/fcr + /cC 



lg(^'^(fcr + fcC) 
2 kr + kC 



-H[^\kT + kO 



(66) 
(67) 
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Now we are ready to consider the creation of the gravitons during the transition from some initial to final states. The 
initial vacuum state |Oin) is determined by ai^jOin) — 0, where Ok is the initial annihilation operator for r^. The relation 
between annihilation and creation operators for the initial and final states is given by the Bogoliubov transformation 



6k = B+ik)a^ + B_{kya\ , 
S+(fc)*4 + 5-(fc)a-k 



(68) 
(69) 



where l-B+P — = 1. Because we are working in the Heisenberg description the vacuum state does not change 

during the evolution. It results that 6k|0in) = B^{k)*a^_^^\Oin) is differ from zero when B^{k)* is a nonzero function. 
This means that in the final state graviton field considered is no more in the vacuum state without particles. The 
number of produced particles in the final state is given by 



|Oin> =|i?-(fc)|' 



Using relations ([44]) and (|45)) and the Bogoliubov transformation ((68)) and (|69|) we obtain 
Mk,Ti)gf{k,Tf) - gi{k,n)ff{k,Tf) 



B-{k) 
B+{k) 



f}{k,Tf)gf{k,Tf) - g*f{k,Tf)ff{k,Tf) 
Mk,n)gj{k,Tf) ~ gi{k,n)f^{k,Tf) 
ff{k,Tf)g}ik,Tf) - gf{k,Tf)f*{k,Tf) 



i [fi{k,Ti)gf{k,Tf) - g,{k,T,)ff{k,Tf)] , 
-■i [Mk,Ti)g*j:{k,Tf) - g,{k,n)f}{k,Tf)] 



(70) 

(71) 
(72) 



where simplifications come from the Wronskian condition In the calculations we set = ri = —20 and 

Tf = T2 = —1. These boundaries fully cover the region of gravitational waves creation. The energy density of 
gravitons is given by 



dPg 



2-fiw ■ 



\B-ikW 



(27rc)3 

where we used definition (|70p . The expression for the parameter flg^ defined by (|25p takes now the form 



k — V ■ 27raf ■ 



^today 
Of 



where 



he 167r2 

Pc 



167r2- 197.3- 10-15 [MeV-m] 
34 . I032[ni4/s4]l.05 • 10-5 • /i2[GeV/cm3] 



3.66 • /in • 10 



-49 



[Hz 



(73) 



(74) 



(75) 



In the calculations we set present value of the Hubble factor for /lo = 0.7. In Fig. [3] wc show spectrum calculated 
with formula (j74[) . The obtained spectrum is extremely weak in the present epoch. The reason of this tiny amount of 
the background gravitons is the presence of the standard inflationary phase. The super-inflationary phase is placed 
before the inflation so the energy of gravitons decreases about 10^^ times during this further phase. To see better how 
presence on the inflation affect this spectrum we show in Fig. [3] the spectrum of relic gravitons in the model without 
the inflationary phase. 

It is clear that in such a model amount of relic gravitons would be extremely large. As we mentioned before it is 
possible that GUT energy scales cover partially with the inflation. In this situation the present amount of the relic 
gravitons would be higher than this in Fig. [3] and reaches to higher energies. For the present state we do not know 
the duration of the inflationary phase exactly. To estimate this value it is necessary to measure the spectrum of both 
scalar and tensor parts of primordial fluctuations produced during the standard inflation. For the present day we 
know only the contribution from the scalar part. The further generation of CMB telescopes in needful to improve 
the knowledge the properties of inflation. From our calculation we can see however that the inflationary phase is 
necessary. Without the inflation after the super-inflation the present amount of gravitons would be easily in reach of 
present observational skills. From this point of view we have found the next motivation to support the inflationary 
model. 



IV. SUMMARY 



In summary, we have calculated the spectrum of gravitons produced during the super-inflationary phase induced 
by Loop Quantum Gravity effects. In the calculations we considered inverse-volume corrections to the dynamics and 
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FIG. 3; Spectrum of relic gravitons for the model with j = 100 and 1 = 3/4. Frequency scale in Hertz. 
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FIG. 4: Spectrum of relic gravitons for the model without inflation. Frequency scale in Hertz. 

to the equation for the tensor modes. We have solved analytically equation for the tensor modes in the quantum and 
classical regimes. Both solutions we had matched by numerical solution for background dynamics. We have obtained 
spectra of relic gravitons for the models considered. In the first model we assumed the presence of the inflationary 
phase after the super-inflation. In that case we obtained presently a negligible amount of relic gravitons. However 
in the second model without inflation the present amount of graviton background would be unnaturally high. This 
results state that the period of inflation after the super-inflation is necessary to avoid the problem of relic gravitons. 
Nowadays the main motivation to introduce the inflationary phase comes from ability to creation of fluctuations. Our 
investigations based on loop quantum cosmology support the inflationary model. 

Results obtained in this work differ from our previous investigations ^ . The difference comes mainly from different 
assumptions about the evolution of the Universe after the super-inflation. The previous results correspond to the 
situation when inflation and GUT energy scales cover each other rather than to the case of very short inflation. In 
this paper we have considered both models agreeing with the present paradigm and its modification. 
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